function [abserr relerr] = cmpder(u, d1u, d2u, params)
  Nx = params.Nx;
  Ny = params.Ny;
  Nt = Nx.*Ny;

  u = u(1:Nt*3);
  d1u = d1u(1:Nt*3);
  d2u = d2u(1:Nt*3);

  u = reshape(u, Nt, 3);
  d1u = reshape(d1u, Nt, 3);
  d2u = reshape(d2u, Nt, 3);

  %std(u)
  %std(d1u)
  %std(d2u)
  abserr = std(d1u-d2u);
  % max(abs(d1u-d2u))
  relerr = std(d1u-d2u) ./ std(d1u);

end
